To appear in The Astronomical Journal 



The Redshift Evolution of the Lyo; Forest 

Tae-Sun Kim^, Esther M. Hu^, Lennox L. Cowie^, and Antoinette Songaila^ 
Institute for Astronomy, University of Hawaii, 2680 Woodlawn Dr., Honolulu, HI 96822 
tsk@ifa.hawaii.edu, hu@ifa.hawaii.edu, cowie@ifa.hawaii.edu, acowie@ifa.hawaii.edu 

ABSTRACT 

We have analyzed the properties of low column density Lya forest clouds (12.8 
< logiVni < 16.0) toward 5 QSOs at different redshifts, using high signal-to-noise 
data from the HIRES spectrograph on the Keck I 10-m telescope. The results 
are used to study the evolution of these clouds in the interval 2.1 < 2; < 3.5. 
We find: 1) The differential density distribution of forest clouds, /(A'hi), fits 
an empirical power law with a typical slope of —1.35 — > —1.55 for clouds with 
A'^Hi < 10^^'^ cm~^ and changes extremely slowly with redshift over this A^hi range. 
2) The deviation of /(A^hi) from this power-law distribution at higher column 
densities depends on redshift: for higher redshifts, the departure from a power- 
law distribution occurs at higher A^hi? which means that higher column density 
clouds rapidly disappear with decreasing redshift. This may be consistent with 
an invariant shape for the different column density distributions, but with a shift 
to a lower column density normalization for systems at different redshifts driven 
by the overall expansion of the universe. 3) The line widths of the forest clouds 
increase as redshift decreases. 4) The correlation strength of the forest clouds 
seems to increase as redshift decreases. We discuss these results in terms of the 
evolution of the IGM comparing the results to models and analytic descriptions 
of the evolution of structure in the gas. 
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1. Introduction 

Lya forest clouds imprinted in the spectra of high-redshift QSOs provide a unique 
opportunity to study the evolution of cosmic structure and the distribution of matter in the 
universe up to 2; ~ 5. The low column density end of the forest is of particular interest, 
since the number density of forest clouds increases as A^hi decreases down to the current 
observational limit, A^hi = 3 x 10^^ cm"^ (Hu et al. 1995, Lu et al. 1996, Kirkman & Tytler 
1997), and the internal density in these clouds is thought to lie close to the mean baryon 
density in the intergalactic medium (IGM). Although much observational progress has been 
made since the advent of the Keck I 10m telescope (which has provided much more sensitive 
observations of the high- 2; clouds) and HST (which has accessed the low-^ clouds), the 
evolution of the low column density clouds with redshift has not yet been described. The 
reason for this is the need for time-consuming high-resolution spectroscopic observations 
which are essential to detect the forest clouds with Hi column densities less than 10^^ cm~^. 
It is this topic which we address in the present paper. 

Recently, a number of CDM-based numerical simulations have shown that Lya forest 
clouds are a natural consequence of structure formation in the intergalactic gas (Cen et al. 
1994, Zhang et al. 1995, Hernquist et al. 1996, Dave et al. 1997, Zhang et al. 1997). Ahhough 
the detailed simulation results and input cosmological parameters vary from model to model, 
there are general trends which may be understood in broad analytic terms {e.g.., Hui et al. 
1996, Bi & Davidsen 1997). In particular, the low column density forest (A^hi < 10^^ cm~^) 
forms in regions where the deviations from the mean density are still small, and shocking 
has yet to occur. (The lowest column density clouds may arise in regions with density 
lower than the mean density.) The normalization of the column density distribution in this 
range depends on the ionizing flux, the mean baryon density, and the ionization history 
of the gas, but the slope of this distribution depends primarily on the amount of power 
in the fluctuation spectrum at near-Mpc scales, and weakly on the thermal history. Thus 
the evolution of the distribution function in column density of these clouds is a powerful 
diagnostic of the formation of structure in the IGM. 

Here we use high-resolution, high signal-to-noise spectroscopic data obtained from the 
Keck I 10m telescope to study the physical properties of lower column density (A^hi < 
10^^ cm~^) forest clouds at z > 2 in three different redshift bins. In particular, we ad- 
dress the following three questions: 1) Is there any change in the number density of the 
lowest column density forest clouds with redshift? 2) Is there any change in the broadening 
of the clouds with redshift? 3) Is there any significant change in clustering properties with 
redshift? Sec. ^ describes the observations and data analysis. The observational results are 
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given in Sec. 0. In Sec. ^ we discuss the observations from the perspective of the numerical 
simulations and also in a broader context. We summarize our conclusions in Sec. ^. 

2. Observations and Data Analysis 

The spectra of five QSOs (Q1422+231, Q0014+813, Q0302-003, Q1700+643, and Q1623+268) 
suitable for studying the Lya forest at z > 2 were selected from an on-going program to 
study QSO absorption line systems (Songaila 1997). The properties of the clouds seen toward 
Q0014+813 and Q0302-003 have already been described in Hu et al. (1995), but the cur- 
rent analysis includes subsequently obtained data which substantially increased the exposure 
times. All QSOs were observed at the Keck I 10m telescope with the HIRES spectrograph 
at a resolution, R, of 36,000 (c/. Cowie et al. 1995, Hu et al. 1995). The properties of the 
quasars and the exposure times are summarized in Table 0. 

The wavelength ranges used for this study were chosen to exclude any higher order 
hydrogen Lyman absorption lines beyond Lya and to avoid the proximity effect. Two of the 
QSOs, Q1422+231 and Q0014-I-813, have partial Lyman limit systems within the selected 
wavelength ranges. Lya forest clouds within ±1000 km s~^ of the partial Lyman limit 
systems in both QSOs were excluded in order to sample only the "quiet" forest (without 
possible influence from the partial Lyman limit systems). Table ^ also gives details of the 
wavelength range and cloud properties. 

Voigt profiles, which at these low H I column densities are Gaussians, were fitted to the 
absorption lines using an automatic line fitting program (Hu et al. 1995). Because the lines 
are formed by a combination of kinematic and thermal broadening, the profile fitting method 
should not be viewed as a direct model of the underlying reality, but rather as an extremely 
powerful parameterization of the spectrum which can be compared to similar analyses of 
simulated model spectra {e.g., Dave et al. 1997, Zhang et al. 1997) to test the models. The 
present fitting program selects lines with a central optical depth, r > 0.05, corresponding to 
an Hi column density, Ahi = 2 x 10^^ cm~^ at 6 = 30 km s~^ (for thermal broadening the 
Doppler parameter, b = ^2kT/mp, where T is the temperature for a thermally broadened 
line, rUp is the proton mass, and k is the Boltzmann constant), and then fits the spectrum 
with an ensemble of profiles characterized by the Hi column density, the line width b, and 
the central redshift z. The program determines a minimal ensemble of lines, since lines can 
always be fit by a larger number of narrower components. This ambiguity can result in 
significant differences and may be the primary cause of the differences between Kirkman & 
Tytler (1997) and Hu et al. (1995), who find somewhat different b value distributions at the 
same redshift. The present analysis is self-consistent in that the five quasars are examined 
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using the same techniques, and therefore the measured evolution should be robust. However, 
we discuss in Sec. p.5| alternative techniques which show the basic results in other forms. Our 
evolution also matches well to the results of Lu et al. (1996) at higher z, despite the possible 
differences in analysis techniques. 

All QSO spectra, except for that of Q1623+268, have S/N of ^40 per resolution element, 
and Q1422+231 has S/N of near 140 per resolution element. In the case of Q1623+268 the 
low blue sensitivity of the HIRES CCD at wavelengths shorter than 4000A dropped the S /N 
to ~ 20 over the wavelengths studied. To construct a homogeneous data set, we only used 
the lines with Nm > 6.3 x 10^^ cm~^ from the 4 other QSO spectra. This also restricts us to 
the column density range where any incompleteness corrections are small (Hu et al. 1995). 
Although it is hard to quantify the saturation of the lines because of the spread in b values, 
we begin to see these effects at about A^hi ^ 1-5 x 10^^ cm~^ for 6 = 30 km s^^ if we define 
a saturated line as a line with a residual flux = exp(— r^,) less than 0.025. As long as the 
lines did not show any severe departures from a Voigt profile, all saturated lines were fitted 
with a single component. 

All possible metal-line candidates were examined. Unlike the case for the z < 1 forest, 
it is not straightforward to distinguish metal lines in the forest at z > 2 due to the higher 
number density of forest clouds. Keeping in mind that most metal lines have lower b values 
(< 15 km s~^) than those of forest clouds when fitted as Hi lines (Hu et al. 1995), it is 
rather easy to recognize Civ (AA1548, 1550) and Si IV (AA1393, 1402). With the exception of 
the spectrum of Q1700-I-643, most CiV and Si IV doublets commonly associated with higher 
column density forest clouds either lie outside the wavelength regions under consideration or 
were not found. Metal lines from the lower redshift metal line systems were also examined 
and excluded. However, broad metal lines such as Si ill (A1206) are hard to identify correctly 
in the forest and may not be completely removed. In any case, the fraction of metal lines 
in the forest will not be higher than 5%, and the number counts are dominated by Poisson 
noise. 

For Q1700+643, which is known to have at least 10 low-redshift metal line systems 
(Reimers & Vogel 1993, Vogel & Reimers 1995, Rodriguez-Pascual et al. 1995), a thorough 
search was made for metal lines and these were removed. However, there are 9 unidentified 
lines with 6 < 20 km s~^ in the Q1700+643 forest compared to 1 line in the Q1623+269 forest 
over the column density range, A^hi = 10^^'^ — 10^^ cm~^. We included these unidentified 
lines in our list of Lya forest clouds; however, they are most likely unidentified metal lines 
from the metal line systems towards Q1700+643. Since the number of unidentified narrow 
lines is at most 11% of the 6 < 20 km s~^ systems, the inclusion or exclusion of these lines 
does not affect the conclusions reached in this analysis. 
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We grouped the 5 QSOs into three different redshift bins (Table |l|), with central values 
of < 2 >= 2.31 (Q1623+268, Q1700+643), < z >= 2.85 (Q0014+813, Q0302-003), and 
<z>= 3.35 (Q1422+231). 

3. Results 

3.1. Differential Density Distribution Function 

The differential density distribution function (hereafter DDF or /(A^hi)) is defined as the 
number of absorbing systems per unit redshift path per unit column density as a function 
of A^Hi- The redshift path X{z) is defined by X{z) = |[(1 + — 1] for go = or by 
X{z) = l[{l + zf/^ - 1] for go = 0.5. 

Fig. |l| shows the DDF, without any correction for blending or incompleteness, at different 
redshifts. The stars represent < z >= 3.35, the triangles < z >= 2.85, and the diamonds 
< z >= 2.31. The open circles are for < z >= 3.70 from Lu et al. 1996. The solid line 
represents the single power- law fit, f{Nm) = 4.9 x WN^^I'^^, taken from Hu et al. (1995) for 
their data at 2; ~ 2.85. The dotted line represents the single power-law fit for z ~ 3.7 Lya 
forest clouds, /(iVni) oc A^h/'^^ estimated from Fig. 6 of Lu et al. (1996), and matched to the 
same /(iVni) at log A^hi = 14.15. This slope of -1.55 at < z >= 3.7 was obtained from the 
fit to the column density range iVni = lO^^-^^^^ '^ cm^^ with the incompleteness correction 
applied. However, over the column density range AThi = 1012-3-14.5 ^^ggQ-^ 
found a slope of -1.46, which matches the value found by Hu et al. (1995) and Kirkman & 
Tytler (1997) at redshifts near 3. 

There are two notable features of Fig. [^. First, the four DDFs roughly follow a single 
power- law, /(iVni) oc N^f, with (3 ~ 1.4 over the column density range Nm = 10 — 
]^gi4.3 ^he DDFs do not show any sign of flattening in this range with the possible 

exception of the < z >= 3.70 DDF at lower Hi column density. However, the lowest 
A^Hi points for Lu et al. have higher levels of incompleteness correction and this result is 
correspondingly uncertain. Earlier results claimed that there was a flattening in the DDF 
at A'^Hi < 10^^'^ cm~2 (Petitjean et al. 1993). However, this flattening is mainly caused by 
blending in low resolution data and was not confirmed from the Keck HIRES data down 
to A^Hi ~ 10^2-5 at 2; ~ 2.85 (Hu et al. 1995). We now see this result holds over the 

2 = 2 — 4 redshift range, and also provides a good description of the Lu et al. (1996) results. 

Second, at the higher column density end, A^hi > 10^^'^ cm"^^ the DDFs for the redshift 
intervals differ substantially. Except at the highest redshift there is a strong lack of forest 
clouds at high column densities compared with the number expected from the power-law 
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distribution at A^hi < 10^^-^ cm~^ (see Petitjean et al. 1993). This deficiency cannot be a 
result of misestimates in tlie column density owing to saturation effects since the total number 
of clouds is deficient relative to the extrapolated power law (Hu et al. 1995). The < z >= 3.70 
data of Lu et al. (1996) show only a slight steepening at the higher column densities, but 
by < 2; >= 3.35 forest clouds deviate from the power-law at A^hi ~ 10^^'^ cm~^, while 
those at < z >= 2.85 and < z >= 2.31 deviate from the power-law at A'^hi ~ lO^'''^ cm~^. 
Furthermore, the degree of deficiency in cloud numbers compared to predictions from the 
single power-law fit also depends on redshift. The < z >— 2.31 clouds deviate from the 
power-law more rapidly than the < z >— 2.85 clouds and the < z >— 3.35 clouds. There 
are no clouds with A^hi > 10^^'^ cm~^ at < z >= 2.31. The redshift path for < z >= 2.31 is 
1.5 times larger than that for < z >= 3.35, so this deficiency of higher column density clouds 
in the lower redshift range cannot be an observational selection effect. Thus the break in 
the density distribution function appears to be strengthening and migrating to lower column 
densities as the redshift decreases. 



3.2. Number density of Lya forest clouds with redshift 

The evolution in Lya forest clouds with redshift has traditionally been explored by 
looking at the number density of clouds in a given column density range as a function of 
redshift. The number density per unit redshift per line of sight is expressed as N{z) = 
No{l + z^, where A'o = ^oT^R^cH^^ . The parameter $0 is the local comoving number 
density of forest clouds, -kR^ is the local cross section of a cloud, and c and are the speed 
of light and the Hubble constant, respectively. This yields 7 = 1 for g'o = and 7 = 0.5 for 
qo — 0.5 for non-evolving clouds. However, it is clear from the evolution of the shape in the 
DDF that the measured 7 will be a function of the chosen column density interval and of 
redshift, except at low column density. This greatly complicates this type of analysis, and 
comparisons between individual studies based on this method are difficult. 

HST observations show a value for 7 = 0.48 ± 0.62 at 2; < 1, which is consistent with no 
evolution (Bahcall et al. 1993). However at higher redshift, z > 2, a rather wide range of 7 
values has been obtained, from 7 = 2.75 ± 0.29 (Lu et al. 1991; indicating strong evolution) 
to 7 = 2.36 ± 0.40 (Bajtlik et al. 1988) to 7 = 1.89 ± 0.28 (Bcchtold 1994; indicating weak 
evolution). An inflection in N{z) from the rcdshift-limited samples has also been pointed 
out, suggesting that 7 is different in different redshift ranges (Lu et al. 1991, Bechtold 1994, 
Impey et al. 1996). These results have been used to suggest that Lya forest clouds do not 
consist of a single population and that a different population becomes dominant at different 
redshifts (Boksenberg 1995). 
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Keeping in mind that N{z) is sensitive to the adopted threshold A^hi (a reflection of the 
form of departure from the single-power law distribution at A^hi > 10^^'^ cm~^), we fitted the 
data to a conventional single power-law form, for lines with A^hi = 10^^'^^ — 10^^ cm~^. This 
corresponds roughly to the conventional equivalent width threshold, W = 0.32A, and enables 
us to compare the present data directly with the results of Bechtold (1994) and Impey et al. 
(1996). Fig. D shows N{z) as a function of redshift. The vertical bars indicate the la errors 
and the horizontal bars indicate the bin sizes for the individual QSOs. No incompleteness 
correction was made; over this column density range, it is rather small (less than 10%) and 
within the la error bars. The filled circles show the Keck HIRES data for our 5 QSOs and the 
open circle that of Lu et al. (1996). The open triangles show the Keck HIRES data only for 
the lower column density systems with A^hi = 10^^'^ — lO^'* '^ cm~^. The diamonds are taken 
from Impey et al. (1996) for HST observations of PG1222+228 and PG1634+706 and the 
open squares are taken from Bahcall et al. (1993) for HST observations of PKS 0044+03, 3C 
273, and 3C 351. The region shown by dashed lines represents the measurement of Bechtold 
(1994), which is a factor of 2 lower than our observations. This is probably explained by 
the 7 times higher resolution in our data, which is capable of deblending neighboring lines 
which, especially at higher redshifts, can be a serious problem. Our data yield a factor of 
2.3 more lines at z ~ 3.2 and a factor of 1.3 more lines at z ~ 2.4 than found by Bechtold 
(1994), consistent with this interpretation of the discrepancy. 

The maximum likelihood method (Murdoch et al. 1986) was used to derive the best 
estimated 7 values for the above data. The best fit using only the present Keck data in our 
study at 3.5 > 2 > 2 (the dotted line) yields 7 = 2.78 ±0.71, consistent with 7 = 1.89 ±0.28 
by Bechtold (1994) within a 2a error. The best fit for A^o is A"o = 3.12. The uncertainty 
in A'o was not calculated since it is strongly correlated with 7. At lower redshifts, z < 1.5, 
line blending becomes less of a problem than for z > 2 systems, since the number density of 
clouds is smaller, and the HST counts should be directly relatable to the Keck counts. When 
our data is combined with the HST observations (Bahcall et al. 1993, Impey et al. 1996), 
the best fit for < 2 < 3.6 yields 7 = 2.15 ± 0.21 and A'^o = 6.89 (the sohd line). Because 
there are only 15 forest clouds at 2; < 1, the maximum likelihood method puts high emphasis 
on the numerous and more densely covered (in redshift space) high-redshift observations, 
resulting in a higher 7 than that obtained with a least-squares fit, 7 = 1.74 ± 0.14. The 
7 value driven by the maximum likelihood method is also higher than the one (7 = 0.48) 
estimated from the z <1 data (Bahcall et al. 1993) or the one (7 = 1.26) from the < z < 4 
data obtained by Bechtold (1994). The latter result is a consequence of Bechtold's lower 
number counts. The value of A^o is also correspondingly smaller (by a factor of 3) than that 
obtained by Bahcall et al. (1993) (A^o ~ 19.0) because of the higher 7. The present high 7 
values show that the number of higher column density clouds aX z > 1 declines rapidly with 
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decreasing redshift. When the < z >= 3.7 systems are included, the value of 7 for < z < 4 
is 2.41 ± 0.18. However, considering the different treatment applied to the different data 
sets, this estimated value should be taken with caution. 

When we count only Lya forest clouds in the column density range of A^hi = 10^^'^ — 
10^^ cm~^, where the DDF is well- fitted by a common power-law for systems at all redshifts, 
we obtain a best fit (the long-dashed line) 7 = 1.29 ± 0.45 and = 39.44 for the z > 2 
systems. This 7 value is consistent within the errors with a non-evolving model (7 = 1), as 
we would expect from our discussion of the DDF. Although the lack of lower column density 
forest clouds at 2; < 1 prevents extrapolation of these fitted values, it appears that the DDF 
and N{z) of these lower column density clouds do not evolve with redshift at z > 2. 



3.3. The change in b value with redshift 

The line widths are characterized by the Doppler parameter, b, in the profile fitting. 
Figs. 3a-c show b — logA^ni diagrams at average redshifts < z >= 2.31, < z >= 2.85, and 
< z >= 3.35. As is shown in the diagrams, most lines with 6 < 15 km s~^ are metal lines. 
There seems to be a reasonably well-defined minimum b value (be) over the A^hi range at each 
redshift when the metal lines are discarded (see also Hu et al. 1995 and Lu et al. 1996, though 
this result is questioned by Kirkman & Tytler 1997), implying that the b value distribution 
is truncated at the low end. (Note the increase in the number of lines with A^hi > 10^^ cm~^ 
with increasing redshift.) Also the number of clouds with 6 < 20 km s~^ (the b^ value at 
z ~ 2.85 shown as a dashed line) decreases as the redshift decreases, suggesting that be is 
increasing with decreasing redshift (see also Lu et al. 1996). 

In order to quantify the changes in be and determine the corrections to the DDF for line 
blending and incompleteness and the profile fitting methodology, we followed the procedures 
of Hu et al. (1995). We generated artificial spectra and fitted them in the same way as we 
treat the observed data. We assumed the slope in the DDF to be 1.5 in order to compare to 
the < z >= 2.85 clouds (Hu et al. 1995), and used a random distribution of clouds at a given 
redshift range. For b values, we used a Gaussian distribution, but truncated it at a value lower 
than the cut-off value, b^, at each redshift. The fitted artificial counts are then compared 
to the input counts to determine the changes caused by blending. This incompleteness 
correction at < z >= 2.31 and at < ^ >= 3.35 as a function of column density is given in 
Table |^, while the incompleteness correction at < 2; >= 2.85 is taken from Hu et al. (1995). 
Fig. 4 shows the incompleteness corrected differential density distribution functions. Since 
the line density is larger at higher redshift, the incompleteness correction at < 2; >= 3.35 is 
slightly more sensitive to the assumed P and the fitted Nm range. Assuming P = 1.5, the 
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resulting /3 is ~ 1.58 for A^hi = lO^^'^ - lO^^-^ cm-^ and ~ 1.53 for Nm = lO^^.s _ ^qm.s ^^^-2^ 
The incompleteness correction aX < z >= 2.31 is not as sensitive to the assumed f3 as the 

< z >= 3.35 clouds. Here, the observed low (3 of 1.28 results in a low f3 of 1.35 after the 
incompleteness correction. For < z >= 3.35, f{Nm) = 3.3 x lO^iVH^-^'^^^"-^^ {2a error). For 

< z >= 2.31, f{Nui) = 1.5 X 10''A/'h/-^^^°-°^ {2a error). We note that using only the two 

< z >=2.85 systems, Q0014+813 and Q0302-003, gives an incompleteness corrected DDF 
having a best fit of f{Nni) = 5.7 x 10^ N^^-^^^^-^^ {2a error) over the column density range 
A'hi = 10^^'^ — 10^^'^ cm~^. This is slightly shallower than the —1.46 slope obtained by 
Hu et al. (1995) with two additional quasars included for this redshift interval. Although 
there might be some deviation at the lower iVni end in the distribution function (where the 
incompleteness correction is very uncertain), the three DDFs show an excellent match at 
Nyii = lO^'^'^ — 10^^'^ cm~^ where the incompleteness correction is not high. Within the 2a 
errors and the systematic errors in the profile fitting methodology, these (3 values could be 
treated as consistent with the < z >= 2.85 forest clouds with P = 1.46 ± 0.07, but overall 
there does seem to be evidence that there is a very weak flattening of f3 with decreasing 
redshift. However, we conclude again that unlike the case for high column density clouds 
(-^Hi > 10^^'^ cm~^), the DDF does not change rapidly with redshift, within the observational 
and simulated errors for A^hi < 10^^'^ cm~^. 

The procedure of incompleteness correction also produces an estimated 6c value. With 
an assumed f3, three parameters of a truncated Gaussian b distribution are varied until a 
satisfactory match to the observations is obtained: median(6), la{b), and 6c- Fig- H shows the 
averaged simulated b distribution (the solid lines) with the observed b distribution (the dotted 
line) ai < z >= 2.31 and < z >= 3.35. With the large number of lines at < z >= 3.35, the 
cut-off b value, 6c; is rather straightforward to determine: 6c = 17 km s~^ (or = 1.8 x 10^ 
K, if we assume a thermal broadening), Icr = 8 km s~^, and median(6) = 26 km s^^. No 
perfect fit was found for the < z >= 2.31 forest clouds. With 6c = 22 km s^^, more than 
50% of the simulations generated too many lines with 6 = 19 — 22 km s~^ when compared to 
the observations. Assuming that most narrow lines are metal lines (90% of the unidentified 
lines with 6 < 20 km s~^ are from the Q1700-I-643 with 10 metal line systems), we took the 
best fit value for the observation namely at 6 > 20 km s~^: 6c = 24 km s~^ (or = 3.5 x 10^ 
K, if thermally broadened), la = 12 km s~^, and median(6) = 30 km s~^. 

These results show that there is an increase in the 6c value as redshift decreases which is 
also consistent with 6c = 15 km s~^ determined at z = 3.4 4.0 by Lu et al. (1996). These 
results are shown in Fig. p^. Fig. p5| shows the median 6 values for each QSO. The solid line 



is for Nm = 10^^-^ - 10^^ cm'^ and the dotted line is for A^m = 10^^'^ - 10^^ cm'^. Table || 
lists the median 6 values for each QSO. This rising trend in both 6c and the median 6 values 
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suggests that the internal temperatures and/or the level of kinematic broadening are rising 
with decreasing redshift, as we shall discuss further in Sec. ^. 



3.4. Correlation Function 

The simplest way of testing inhomogeneity in the space distribution of Lya forest clouds 
is to calculate the two-point velocity correlation function, ^. Among different ways of calcu- 
lating ^, the most widely used method is comparison between the observed number of pairs 
(Nobs) and the expected number of pairs (Nexp) from the random distribution in a given 
velocity bin (Av); ^{Av) = Nobs{Av)/Nexp{Av) - 1, where Av = c{z2 - zi)/[l + {z2 + zi)/2] 
and Zi and Z2 are the redshifts of two lines. Although this formula is correct only in the limit 
of an ideal large sample and is sensitive to the bin size, general information on clustering 
can be obtained by this method. 

Searches for clustering in Lya forest clouds have been a function of resolution, S/N, bin 
size, and possibly line of sight, and there have been rather mixed results to date. Sargent 
et al. (1980) and Ranch et al. (1992) find no positive correlation at any velocity scale. On 
the other hand, Cristiani et al. (1995), and Hu et al. (1995) show an anti-correlation at 
Av = 600 - 900 km s~^ as well positive correlation at Av < 200 km s ^. While this 
anti-correlation is also claimed to exist by Meiksin & Bouchet (1995), Cristiani et al. (1997) 
failed to confirm the reality of anti-correlation from a larger data set, but did confirm the 
positive correlation at Av < 300 km s~^. 

The correlation functions at the three redshifts were calculated relative to the average 
of 500 simulations with a 50 km s~^ velocity bin for each QSO. We used lines with A^hi < 
10^^ cm~^ for this calculation, although there are few lines with A^hi > 10^^ cm~^. 



Fig. shows the averaged correlation function as a function of redshift for systems 
with column densities above a threshold value, A'ni^th, such that logA'nLth > 13.8. The 
dotted lines are the Icr Poisson errors. While there is suggestion of a positive correlation 
for < z >= 2.31 at Av < 300 km s"^ (a 2.8cr significance at Av = 75 km s~^), there is no 
significant correlation strength at < z >= 2.85. This is in disagreement with what was found 
by Hu et al. (1995) with an overlapping HIRES data set. We will discuss this point later. 
At < z >= 3.35, a positive correlation may exist at a 2.54cr significance at Av = 75 km s~^. 
Given the large time gap between z = 3.35 and the present day, the interpretation of the 
correlation functions with redshift is ambiguous. However, the clustering should evolve 
strongly with redshift aX z = 2.31 ^ 2; = if Lya forest clouds follow more or less the highly 
structured present day galaxies (Hernquist et al. 1996, but see Morris et al. 1993). Fig. ^ 
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shows the correlation strength at Av = 75 km s~^ as a function of threshold column density 
^Hi,th5 for three redshift bins. As Cristiani et al. (1997) pointed out, aX < z >= 2.31 and 
< z >= 3.35, the correlation strengths seem to increase as A^Hi,th increases. At < z >= 2.85, 
this trend does not hold tightly. However, as Fig. ^ shows, when averaged for all redshifts, 
the correlation strength does increase as A^Hi,th increases. 



3.5. Alternative Descriptions 

Figs. and [SB show two alternative descriptions of the data {e.g., Croft et al 1997). 



Fig. shows the normalized cumulative distribution of the flux decrements D = 1 — e'^ for 
the five quasars computed in the wavelength ranges specified in Table |l]. The distributions 
agree well with similar observations by Ranch et al. (1996) and also with model predictions 
for a suitable choice of Qkh^T-^-^J {e.g., Ranch et al. 1996, Croft et al. 1997). Fig. ^ 
shows the number of times the specrum crosses down across a given decrement per unit 
X. In the absence of blending and given the invariance of the DDF, narrower b values will 
result in a higher number of crossings. Thus, the increased rate of crossings seen at higher 
z for the larger flux decrements is an alternative method of seeing that the lower z clouds 
are broader. At small flux decrements blending becomes important, and because this effect 
is larger at higher z, this reverses the sign of the effect. 



4. Discussion 

4.1. The column density distribution function 

The number density per unit redshift path X per unit Nm (the density distribution 
function) as a function of redshift depends on the history of the UV ionizing background, the 
mean baryon density, and the underlying cosmology, and only detailed numerical simulations 
can give complete descriptions of these combined effects to compare with the observations. 

Simulations by Miicket et al. (1996) show broad consistency with the present data, with 
high column density clouds {Nm > 10^'^ cm~^) evolving as 16.8(1 + z)^^^ + 0.687(1 + z^-^ 
while low column density clouds (A^m < 10^^ cm~^) evolve only as (1 + z)^'"^ . However, Fig. 
^ shows that the observed rise for the high column density (10^"^ < Ahi < 10^^ cm~^) clouds 
is somewhat steeper than this prediction. Lu et a/.'s results aX < z >= 3.7 are also shown, 
and are consistent with the fit obtained for < z >= 2.31 3.35 . The more sophisticated 
treatment of Dave et al. (1997), using Voigt profile fitting to artificial spectra generated from 
their SPH simulations, shows a DDF that drops by a factor of 1.5 — 2 for 2; = 3 — > 2 at 
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-^Hi > 10^^'^ cm~^, while at A^hi < 10^^'^ cm~^, the number of clouds increases by a factor 
of 2 for z = 3 ^ 2. These features in the predicted DDF do not match the observations, 
which show a larger number of clouds at 2; ~ 3 compared to those at z ~ 2 only for those 
clouds with A^Hi > 10^^'^ cm~^, but this may be due at least in part to the profile fitting 
methodology. Both Dave et al. (1997) and Weinberg et al. (1997) find relatively steep profiles 
with (3 = 1.7 as do Miicket et al. (1996). However, Zhang et a/.'s (1997) similar analysis of 
their hierarchical grid code simulation shows a remarkably good agreement both with the 
overall shape of the DDF at z = 2.9 {f3 ~ 1.5), and with the redshift evolution seen in the 
present data, finding 7 for 14^ > 0.32A selected clouds of 2.95 ± 0.4 for z = 1.5 — > 3.5 and 
0.24 ± 0.03 for z < 1.5, which is in good agreement with Fig. |^ and the discussion of Sec. 
|3.3| . It is unclear whether the differences between Dave et al. and Zhang et al. (which are 
also reflected in the b value distributions discussed below) are a consequence of details of the 
codes, choice of parameters (particularly the fluctuation spectra), or details of the profile 
fitting analysis, but the results do suggest that within the tunability and uncertainty in the 
current models that it is possible to obtain remarkably good agreement with the data. 

The analytic models of Hui et al. (1996) give considerable insight into the origin and 
invar iance of the slope. At these low column densities the temperature is set by the balance 
between adiabatic expansion and ionization heating, and may be described by a polytropic 
equation of state T ~ p'^'^, where (7 — 1) lies between and 0.62, depending on the ionization 
history. The slope of the DDF is then given by 



P = l + 



0.96 - 2 ((To - 1) 
1.68-0.7(7-1) 



where cTq is the RMS density fluctuation at scales close to the Jeans length (Hui et al. 
1996). For the observed (3 values and for the possible range of 7 values, cTq is extremely 
tightly constrained to a value of 1.15, with a plausible range of only ±0.1. The flattening 
of the P value with decreasing redshift would require 7 to be closer to isothermal at the 
lower redshifts, which is counter to the expected evolutionary trend. This may, however, be 
related to reionization of the He if this occurs around z = 3 as has been suggested based on 
the evolution of metal line ratios (Songaila & Cowie 1996). 

The deviation from the power law at the high column density end would then correspond 
to the overdensity at which shocking occurs {e.g., Miicket et al. 1997), which is approximately 
ten times the mean density at the given redshift. If the column density at which the deviation 
occurs does correspond to the shocking overdensity then we expect it to scale roughly as 
(1 + z)^ T^^'"^ J^^ , where T is the temperature and J the ionizing flux. If the ionizing flux 
is roughly constant over the redshift range (Haardt & Madau 1996) and the temperature 
changes only slowly with redshift as in the limit of early ionization (Hui & Gnedin 1996), A^dov 
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would drop by 0.8 in the log between 2; = 3.5 and z — 2.3, which is somewhat larger than the 
observed drop from log A^hi — 14.8 to 14.3. This suggests that J is falling with decreasing 
redshift or that the temperature is rising. The latter effect could again be associated with 
helium reionization. 



4.2. The evolution of h values 

If the LycK forest lines were intergalactic clouds in thermal equilibrium with the meta- 
galactic UV ionizing background without any other energy source, the line width of the clouds 
(or the Doppler parameter, 6) would measure the temperature of the clouds. However, within 
the structure evolution models the low column density forest clouds (A^hi < 10^"^ cm~^) are 
primarily kinematically rather than thermally broadened. Two important observations may 
be taken from the Keck HIRES data. First, there seems to be a rough lower cut-off h value 
(6c) which increases with decreasing redshift. Most forest clouds have a h value similar to or 
greater than 6c, whereas most lines with b < be are identified as metal lines. Hu et al. (1995) 
interpreted 6c as the minimum 6 value in a single region, with most Lyo; clouds being a blend 
of several components. Secondly, taking the median 6 as 34 km s~^ and 6c = 20 km s~^ sug- 
gests that the additional broadening due to the multiple components is very similar to h^. 
(Cowie et al. 1995). The increase in the h values with decreasing redshift is likely to be a 
reflection of the ongoing increase in structure. 

The numerical simulations of Dave et al. (1997) give higher 6 values than are observed 
at ^ ~ 3, while the simulations of Zhang et al. (1997) provide fairly good agreement with 
the observations. (Both are broader than the distribution of Kirkman & Tytler 1997.) This 
is a similar effect to that seen in the DDF, and may indicate some fundamental difference 
between the SPH and grid simulations. There is as yet no discussion of the evolution of the 
6 value distribution, though the overall increase in the width with redshift is in the sense 
expected from the models. 

For the higher column density clouds which have entered the shocked regime the widths 
will reflect primarily the thermal broadening, and the minimum widths may then correspond 
to minimum temperatures [e.g.., Haehnelt et al. 1996). However, even for these clouds the 
evolution of internal temperatures governed solely by ionization equilibrium is insufficient to 
provide an explanation of the evolution of 6c. Songaila & Cowie (1996) show that Siiv/Civ 
in the forest clouds at 2.1 < 2; < 3.8 is well described by a hypothetical Hell ionization front 
at ^ ~ 3.1, with Si/C = 3 times the solar value, since the ratio of Silv/ClV is sensitive 
at the high energy end in the assumed mctagalactic UV flux. If there is a Hell ionization 
front at z ~ 3.1, this ionization causes a change in the mctagalactic UV flux, and thus the 
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b value in the forest clouds. We tested this effect using a broken power law for the UV flux, 
following the same procedure as Giallongo & Petitjean (1994), where the amount of Hell 
ionization is characterized by a parameter B, which is the relative flux at 4 ryd and 4.5 
ryd. Assuming the flux at the Lyman Limit (1 ryd) of Jo = 10~^^ erg cm~^ s~^ Hz~^ and a 
single power-law for the UV energy spectra with a = 1.5 (Jy oc i^~") (Cooke et al. 1997), we 
ran the photoionization model of Ferland (1993). The calculation was made for clouds with 
A^Hi = 10^^'^ cm~^ and we used a 10~^ solar metallicity with Si/C = 3. Fig. ^ shows the 
model calculations with different ionization parameters, F ~ Jo/'^e- The long-dashed line is 
for B = 1 (no Hell ionization); the dotted line, the dashed line, the dot-dashed line, and 
the dot-dot-dashed line are for B = 10, B = 100, B = 1000, and B = 10000, respectively. 
Songaila & Cowie (1996) found a satisfactory result for the Si and C ratios with B > 100 
for z > 3.1 and B = 1 — 10 for z < 3.1. At < 2; >= 2.85, the model-calculated b for the 
hydrogen clouds is well-matched with the 5 = 1 — 10 spectra at F = 10~^ — 10~^'^ and at at 
< z >= 3.35 and 3.7, the be values are also in the regime of the photoionization prediction 
with B = 100. However, at < z >= 2.31, the b^ value is too high for the photoionization 
model to obtain (for either be = 22 km s~^ or be = 24 km s"^). Thus evolution in b^ 
cannot be fully accounted for by evolution of an equilibrium photoionization temperature. 
A hypothetical Hell ionization at 2; ~ 3.1 can explain the evolution of be values a.t z > 2.5, 
but fails to explain the 6c of 24 km s~^ at < z >= 2.31. The lower redshift clouds clearly 
require a higher dynamical energy input than the high redshift clouds. This may reflect a 
higher level of overall turbulence at the lower redshifts, as evidenced by the increase in the 
median b values. 



4.3. Lya forest contribution to fib 

One of the most important cosmological parameters is the ratio of the baryon density to 
the critical density in the universe, fib, which can be directly compared with the predictions 
of big bang nucleosynthesis. At 2; = 2 — 4 the intergalactic gas may still comprise most of 
the baryons in the universe and thus the QSO absorption line systems can provide a direct 
measurement of fib. 

Numerical simulations of the Lyman forest have generally favored a high fib/i^ (^0.0125) 
{e.g., Ranch et al. 1996, Weinberg et al. 1997), though a more direct study of the Hell 
opacity suggests a lower fib (Hogan et al. 1997). All of these methods are dependent on 
assumptions about the evolution of the shape and normalization of the ionizing field (as well 
as its fluctuations), which determine the cloud internal temperatures and the conversion 
from neutral to total hydrogen. It is also possible that the current simulations may not be 
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adequately dealing with the structure on the smaller scales (Croft et al. 1997), which again 
could result in uncertainties in the relative normalization of the hydrogen optical depth 
relative to f2b. Here we shall adopt a more empirical approach to determining the evolution 
of fib from the observations. 

The contribution to fib(Hi) from the Lya forest can be directly obtained by integrating 
through the DDF. In units of the local critical density pc-, the mass density of Hi in the Lya 
forest is 

^],(HI) = T""^" Nnif{N^i)dNni, (2) 

{e.g., Lanzetta et al. 1991). The parameter mn is the hydrogen mass, c is the speed of 
light. Ho is the Hubble constant, and A^min — A^max is the Hi range of the forest clouds. The 
parameter p is the mean mass per proton mass and is taken as 1.3. 

Because the DDF is everywhere shallower than a -2 power-law slope, fib (Hi) is always 
dominated by the highest column density clouds. In Table ^ we summarize fib (Hi) from 
clouds below 2 x 10^^ cm^^, which is always below the break, and also the total fib(Hl) 
for A^Hi < 4 X 10^^ cm^^. In the lower Nm range which is below the break, fib(Hl) is 
roughly constant at 6 x 10~^ over the redshift range, reflecting the rough invariance of 
the DDF, while in the higher column densities fib(Hl) drops from 2.5 x at z = 3.7 

to 1.1 X 10"'' /i"^ at z = 2.31, reflecting the decreasing number of clouds. 

In order to translate fib(H l) to fib we need to know the appropriately weighted ionization 
fraction F = Hl/H in the clouds. At the high column density end (A^hi-^IO^^ cm~^) the 
ionization state of the clouds is directly probed by metallic line ratios and Cll/Civ ratios, 
around 0.01 (Songaila & Cowie 1996), give typical values of the ionization parameter U = 
0.002 0.03, where U = 4.2 x 10^^ J2i/n in terms of n, the total hydrogen density, and 
J21, the ionizing flux in units of 10~^^ ergs cm~^ s~^ Hz~^ sr~^. Here we have assumed 
Ji, ~ with a = +1.5, following the discussion of Donahue & Shull (1991). The exact U 
value determined from the Cll/Civ ratio depends on the details of the cloud temperature 
(Haehnelt et al. 1996), and the range of U parameters reflects this uncertainty as well as 
the observed range in Cll/Civ ratios. For pure photoionization equilibrium F = 3.4 x 
10~^ f/-i-066 j^Yie corresponding F lies in the range 1.4 x 10~^ —>■ 2.5 x 10~^. As Haehnelt 
et al. (1996) have emphasized, the temperature spread caused by shocking in these clouds 
may result in additional collisional ionization contributions, which can result in the hydrogen 
fraction's being substantially underestimated by photoionization equilibrium. However, the 
b values provide an upper limit to the temperatures of the clouds since 

T < 10^ f (3) 

V13kms"V 



- 16 - 



and we can see from the distribution of b values that there are relatively few high tempera- 
ture clouds in the sample. (At z = 3 less than 19% of the clouds can be above 10^ K.) We 
have integrated the coUisional ionization equilibrium though the cloud temperature distri- 
bution using the assumption that the clouds are thermally broadened to obtain a mean ratio 
A^(H)/A^(Hl) = 1.6 X 10^ at z = 3. It should be emphasized that this is an extreme upper 
limit, since it ignores blending and kinematic broadening contributions to the b distribution. 
The value is smaller at higher redshifts where the b values are lower. 

We have no direct measure of the ionization parameter for the lower column density 
clouds, and instead must resort to indirect arguments. In the linear regime, provided reion- 



ization took place early, T 



n 



+0.6 



and then in ionization equilibrium thi 



n 



+1.1 



(Hui 



et al. 1996, Croft et al. 1997) which, since the b value - density distribution is also relatively 
invariant as a function of nni, implies A^hi ~ n^^'^. The Hi fraction as a function of A^hi 
column density then scales as 



J'. 



AT, 



HI 



21 



In combination with the Hi power law indices of 1.38 — > 1.55 this then implies that the 
contribution to fib is very weakly dominated by the higher column density clouds just below 
the break in the DDF, whence 



C pc 
c Pc 



HI 



/devAT, 



dev 



:i.6 - /3)Fdev 



1 - 



dev 



dNui 

1.6-/3' 



(5a) 
(5b) 



where /dev is the value of the DDF and Fdev the value of the neutral hydrogen at the Adev 



column density. For (3 



1.45 and A",, 



fib(HI) 



= 0.01 A'dev "we obtain 
1.2 X 10-^ h-^ / A", 



dov 



2 X 1014 cm-2 



0.5 



- -^dev ^ -^dev 

from Table ^. Raising (3 to 1.55, the constant rises to 2.1, while reducing N„ 
increases it to 2.43, and the combination would increase it to 2.7. 



(6) 

to 0.001 A'dev 



In order to minimize the downwards extrapolation from the metal line measurements at 
Ahi ~ 10^^ cm~^ we restrict ourselves to the highest z range, 
at the column density range probed by the metals. For Adev 
= 1.2 X 10-'^h'\ and (3 = 1.53 (Table 1), and setting N^i^ 
contribution from lower column density clouds is then 



lere Adev ~ lO^"^'^ cm ^ lies 

^Ql4.8 ^j^g^g fib(Hl) 

_ j^Qii.8 Q]2i~^, the baryonic 



fib = 0.003/1" 



dev 



10' 



(7) 
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In order to obtain a conservative upper bound we assume that the C ll/C IV ratios discussed 

above correspond to clouds with a A^(Hl)~ 10"^^ cm^^ and adopt U = 0.03 at this A^(Hl). 
Scahng back to lO^"^'^ then gives F^^v from photoionization of 5 x 10~^. In combination 
with the upper bound estimate from coUisional ionization, we obtain Fdev > 3 x 10~^ and 
Ob < 0.01 /i-^ 

This is just consistent with low end estimates of the lower bound on from Weinberg 
et al. (1997), but the overlap range is surprisingly small. This discrepancy can be directly 
identified with the adopted limits on J by considering the overdensity at A^dev From the 
definition of U 

" = 2xlO-fi;)"V;!§) cm- (8) 



.o.iy V0.5 

so that the overdensity at ionization parameter U is 

n _ 4.5 \V f/ \"VJ2i 



<n> 1^ 0.01 J \l + z) VO.lJ VO.sj' 

Within the models the deviation density corresponds to an overdensity of roughly 10, and 
we may rewrite this as 

so that the present results require J 21 < 0.2, which then relaxes the fl^, h"^ constraints in 
Ranch et al. (1996) and Weinberg et al. 1997 at the expense of an inconsistency with their 
lower limit estimates of J. 



4.4. The correlation functions 

The evolution of the correlation strength with redshift is very ambiguous. As Morris 
et al. (1993) show, the distribution of the local Lya does not follow the local galaxy distri- 
bution closely, but neither is it quite random. However, numerical simulations show that 
Lya forest clouds follow the distribution of dark matter and do show the galaxy distribution 
more closely at 2; ~ 3 than does the observed local distribution of Lya forest clouds. As long 
as the clouds are not randomly distributed, the evolution of their correlation strength with 
time provides a glimpse of structure formation and/or evolution in the early universe. 

While the positive correlation over a small velocity range is marginally significant, an- 
other interesting feature, the anti-correlation, shows contradictory results. The existence of 
an anti-correlation was first claimed for spectra toward Q0055-269 (Cristiani et al. 1995) and 
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was later confirmed by Meiksin & Bouchet (1995). A similar anti-correlation is also shown 
in the correlation function of the Keck HIRES data at < z >= 2.85 (Hu et al. 1995). How- 
ever, with a larger data set, Cristiani et al. (1997) claimed that there is no anti-correlation. 
Interestingly, the correlation function calculated dX < z >= 2.85 for this study (only for 
Q0014+813 and Q0302-003) using the same Keck HIRES data does not show the anti- 
correlation. The results of Hu et al (1995) used 4 QSOs in this redshift interval (the present 
two together with Q0636-I-680 and Q0956-I-120). Most of the signal in the anti-correlation 
arises from Q0636-I-680, and when this QSO is excluded the sample does not show any strong 
indication of the presence of an anti-correlation. 

This result also gives a possible explanation for the contradictory results of Meiksin 
& Bouchet (1995) and Cristiani et al. (1997). Since Meiksin & Bouchet (1995) used only 
individual QSOs to look for an anti-correlation, they were able to find and confirm the early 
result of Cristiani et al. (1995). However, with the larger data set, most peculiar features in 
correlation functions are averaged, and no longer give the anti-correlation shown by smaller 
numbers of QSOs. This issue requires further work with much larger data sets. 



5. Summary 

We have analyzed Lya forest clouds (12.8 < log A^hi < 16.0) towards 5 QSOs in order to 
study the evolution of these clouds at 2.2 < 2; < 3.5, using Keck I HIRES data. We include 
the Lya forest cloud analysis at 3.4 < z < 4.0 given by Lu et al. (1996) from similar quality 
Keck data to summarize the evolution of the Lya forest at 2 < 2; < 4. 

The differential density distribution function (DDF) does not change rapidly with red- 
shift at the lower Hi column density end, being well fitted with power-laws close to ~ 1.4. 
However, at the higher column density end, DDEs depart from a single power-law; at higher 
redshift, the amount of deviation from the power-law is smaller and the deviation column 
density at which DDEs start to deviate from the power-law is larger. The number density 
with redshift including HST observations aX z < 2 indicates mild evolution in the number 
density, which reflects the behaviors of the DDEs. The minimum, mean, and median line 
widths of the Lya clouds tend to increase as the redshift decreases. The correlation strength 
seems to increase as redshift decreases, at a velocity scale less than Av < 300 km s"^. 
The fluctuation power spectrum at near-Mpc scales (c./., Hui et al. 1996) is found to be 
1.15 ±0.1, while the contribution to Jib from forest clouds (12.8 < logA^ni <14.8)at2; = 3.5 
is flb<0.01h-\ 
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Table 1. Fitted QSOs 







-Grange 








QSO V 




(A) 


^Lya 


dX^ of lines 


Q1623+268 16.0 


2.521 5'^ 20™ 


3850 - 4195 


2.17-2.45 


0.939 


68^^ 


Q1700+643 16.1 


2.722 2 40 


3850 - 4195 


2.17 - 2.45 


0.939 


82 


Q0014+813 16.5 


3.387 8 15 


4510 - 4863 


2.71 - 3.00 


1.023 


118 


Q0302-003 17.6 


3.290 10 00 


4510 - 4863 


2.71 - 3.00 


1.119 




Q1422+231 16.5" 


3.620 9 40 


5105 - 5485 


3.20 - 3.51 


1.233 


146 


Note. — Magnitudes and redshifts are taken from Veron-Cetty & Veron (1996). 




^The redshift path is 


defined as X = 0.5[(1 + zf - 1], for Qo = 0. 






''The number of lines 


is for Nm = 10^^-^ - 10^^-^ cm- 


2 






magnitude. 














Table 2. Incompleteness correction^ 








< z>^ 


2.31 


< Z>^ 


3.35 




log A^HI 


Incompleteness 


log/ 


Incompleteness 


log/ 




12.8 - 13.1 


0.68 


-11.37 


0.44 


-11.02 




13.1 - 13.4 


1.00 


-11.68 


0.72 


-11.54 




13.4 - 13.7 


0.92 


-12.04 


0.78 


-11.93 




13.7 - 14.0 


0.89 


-12.58 


0.89 


-12.53 




14.0 - 14.3 


0.90 


-12.94 


0.91 


-12.90 





^For the incompleteness correction for the < z >— 2.85 forest clouds, the 
results of Hu et al. (1995) were used. 
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Table 3. Median b values of the observed QSOs 



QSO 


iVni = 10i3 « - 
(km s' 


10^^ 


Nm = 10^^-^ - 10^^ cm-2 
(km s"^) 


Q1623+268 


41 




36 


Q 1700+643 


35 




28 


Q0014+813 


34 




31 


Q0302-003 


37 




29 


Q1422+231 


30 




27 


QOOOO-003^ 


31 




29 



Note. — Median b values over the indicated column density ranges. 

^From Lu et al. 1996 



Table 4. of neutral hydrogen 



Parameters < z >= 2.31 < z >= 2.85 < z >= 3.35 < z >= 3.70 



-/Vffl.dev (cm 2) 


1014.3 


1014.3 


1014.8 




f{Nm) 


1.5 X lO^A^i/-^^ 


4.9 X lO^A^i/-^^ 


4.5 X lO^AT-ji-s^ 


9.5 X 108 A^h/-^^ 


nb(logiVHi < 14.3) 


5.7 X 10-8 


6.0 X 10-8 


6.3 X 10-8 


7.2 X 10-8 


nb(log7VHi < 15.5) 


1.1 X 10"^ 


1.3 X 10"^ 


2.1 X 10-^ 


2.5 X 10-^ 
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Fig. 1. — Differential density distribution function without incompleteness correction. The 
distribution with H I column density is shown for each of three indicated redshift bins, with 
1(7 error bars, compared with the results of Lu et al. (1996) (open circles). The solid line 
represents the single power-law fit, f{N^i) — 4.9 x 10^ ^m '^^ ■ The dotted hne shows the 
single power- law fit taken from Lu et al. (1996) for < z >^ 3.7, which is f{Nm) oc N^l^^ for 
12.6 < logA^Hi < 16.0. At A^Hi < 10^^ cm~^, the distribution functions at different redshifts 
look very similar, and lack significant flattening. However, at A^hi ^ 10^"^ cm^^, the amount 
of deficiency in forest clouds compared to the expected power-law distribution depends on 
redshift. The deviation from the power law moves to smaller iVni as redshift decreases, which 
means the higher A^hi forest clouds evolve rapidly. 
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□ Bahcalletal. (1993) 

o Impey etal. (1996) 

• This study 

A This study 

(13.1 <Log Nhi< U.Ocm"^) 

r_"j Bechtold (1994) 

o Lu etal. (1996) 



Log Nhi = 13.77 - 16.0 cm" 
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Log (1+z) 
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Fig. 2. — N{z) vs redshift for the present data (filled circles) and the HST observations 
(open squares and open diamonds). The region outlined with dashed lines shows the results 
from Bechtold (1994). The solid line shows the maximum likelihood fit to the HST data 
and the high column density Keck data (slope 7 = 2.15). The dotted fine shows the fit to 
the Keck data only, with a slope of 7 = 2.78. Triangles show the forest clouds over just the 
column density range A^hi = 10^^-^ — 10^^ cm~^, with a long dashed fine fit with 7 = 1.19. 
The la error bars and redshift binning are indicated for each point. 
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Fig. 3a. — b — logA^m diagram ai < z >— 2.17 — 2.45. The dot-dashed hne shows b — 
20 km s~^. Note that there are not many hnes with A^hi ^ 10^'* cm~^. The vertical dotted 
hnes show the adopted minimum A'hi and the value at which the clouds saturate, while the 
dashed line shows the optical depth selection limit. 
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Fig. 3b. — b — logiVni diagram ai < z >= 2.70 — 3.00. The dot-dashed hne shows 6 = 
20 km s~^. Note that there are many more hnes with A'hi > 10^^ cm~^ compared to the 
same range in column density for z — 2.17 — 2.45. 
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Fig. 3c. — b — logiYni diagram a.t < z >= 3.20 — 3.51. The dot-dashed hne shows b = 
20 km s~^. Note the increase in the number of hnes with 6 < 20 km s~^ over this redshift 
range compared to < 2; >= 2.31 and < z >— 2.85. 
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Fig. 4a. — Incompleteness corrected differential density distribution functions with redshift. 
Diamonds are < z >= 2.31 clouds and triangles are < z >= 3.35 systems, la error bars are 
shown for each point. The dotted line is the same power-law fit shown in Fig. |l|. See the 
discussion in Sec. ^]3] for more details. 
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Fig. 4b. — Corrected and uncorrected differential density distribution functions at < z >= 
2.31. Tlie deviation from tlie power-law fit at column densities above iVni = 2 x 10^^ cm~^ 
may be clearly seen. 
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Fig. 4c. — Corrected and uncorrected differential density distribution functions at < 2; >= 
3.35. 
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Fig. 5. — Observed (dotted histogram) and simulated (solid histogram) b distributions. For 
the < z >= 3.35 clouds, the low cut-off b value is ~ 17 km s~^. For the < z >= 2.31 clouds, 
there is no single best-fit assumed distribution for b. Assuming that most narrow fines are 
unidentified metal lines towards Q1700-I-643 and ignoring the lines with b < 22 km s~^, the 
cut-off b of 24 km s"^ which is shown in the figure best fits the data. 
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Fig. 6a. — Evolution of 6c with redshift. 
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Fig. 6b. — Evolution of median b values with redshift. The sohd line is for A^hi = 10^^-^ — 
10^^ cm^^, and the dotted line is for A^hi = 10^^'^ — 10^^ cm~^. There is a general trend of 
increasing b values as redshift decreases for both column density selections, though the effect 
is more pronounced in the higher column density systems. 
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Fig. 7a. — The rcdshift- averaged correlation functions. The dotted Unes show the la Poisson 
uncertainty. At < z >= 2.31 and < z >= 3.35, there is a positive correlation at Av < 
300 km s~\ 
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Fig. 7b. — The correlation strengths at Av — 75 km s ^ with A'ni.th at different redshifts. 
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Fig. 7c. — The redshift-averaged correlation strengths at A^; = 75 km s ^ with -/Vni.th- 
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Fig. 8a. — Normalized cumulative distribution of flux decrements for the 5 quasars. The 
dashed line shows the < z >= 3.35 system, the solid lines the < z >= 2.85 systems, and the 
dotted lines the < z >= 2.31 systems. 




Fig. 8b. — Number of down-crossings per unit X plotted vs. flux decrement. Redshift 
identiflcations as for part a). 
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Fig. 9. — N{z) vs redshift at A^hi = 10^"* cm ^. The filled circles show our data points, 
while the open circle is taken from Lu et al (1996). For Q0014+813 and Q0302-003, the 
N{z) values are nearly identical, and are marked with a single overlapping data point. The 
data show a uniform trend with a steep rise in the number density evolution. The dotted 
line shows a least-squares fit to N{z) oc (1 + z)^-^^. 
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Fig. 10. — Photoionization and h values. The lines were calculated from the photoionization 
model for the assumed Hell ionization at 2; ~ 3.1. See the discussion in Sec. [4.2| for more 
details. 



